Lifshitz-like transition and enhancement of correlations in a rotating bosonic ring lattice 
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We study the effects of rotation on one-dimensional ultra-cold bosons confined to a ring lattice. For com- 
mensurate systems, at a critical value of the rotation frequency, an infinitesimal interatomic interaction energy 
opens a gap in the excitation spectrum, fragments the ground state into a macroscopic superposition of two states 
with different circulation and generates a sudden change in the topology of the momentum distribution. These 
features are reminiscent of the topological changes in the Fermi surface that occurs in the Lifshitz transition in 
fermionic systems. The entangled nature of the ground state induces a strong enhancement of quantum corre- 
lations and decreases the threshold for the Mott insulator transition. In contrast to the commensurate case, the 
incommensurate lattice is rather insensitive to rotation. Our studies demonstrate the utility of noise correlations 
as a tool for identifying new physics in strongly correlated systems. 



Ultra-cold gases loaded in an optical lattices are becom- 
ing one of the most exciting platforms for exploring com- 
plex phases of strongly correlated systems. Flexible exper- 
imental control over the parameters of the lattice and the 
strength of the interatomic interactions have led to the real- 
ization of fermionization of bosons, i.e. the Tonks-Girardeau 
(TG) regime |1, 2], and a one-dimensional Mott insulator 
(MI) state [1]. Recent experimental advances |3] are open- 
ing a new arena for the investigation of persistent currents 
and various novel physics that emerge when cold atoms are 
trapped in ring-shaped optical lattice. In these systems it is 
now possible to create an artificial magnetic field by rotating 
the atomic cloud lHH^HHtl and to study the exotic phenomena 
that emerge when interacting particles are subjected to strong 
gauge fields. 

In this paper, we investigate the effects of rotation on the 
ground state properties of bosonic atoms confined in a one- 
dimensional ring lattice. When the lattice is rotated with an- 
gular frequency Q,, the Coriolis force generates an effective 
vector potential, and hence a net circulation or non zero vor- 
ticity. The physics of the interacting quantum many-body 
system differs considerably between the commensurate case 
(CO: where the number of atoms, N, is an integral multi- 
ple of the number of lattice sites, L) and the incommensu- 
rate case, ICO. Recent studies jgt] have shown that at a critical 
frequency of rotation, fie, the ground state of CO lattices be- 
comes a macroscopic superposition (cat state) of two states of 
opposite circulation jlOll . Here we focus on the role of the 
interatomic interaction. For CO lattices, we show that a max- 
imally entangled state exists at an infinitesimal value of the 
on-site interaction parameter, accompanied by a discontinu- 
ous rearrangement of momentum distribution, the opening of 
a gap in the energy spectrum and non-analytic behavior of the 
compressibility. These effects are analogous to the sudden 
change in the topology of the Fermi surface in the conven- 
tional Lifshitz transition in fermionic systems ifTTll . in which 



a gap and a discontinuous change in the momentum distribu- 
tion are induced as the pressure approaches a critical value. 
We show that the entangled nature of the ground state as the 
system approaches ilc is accompanied by a large depletion of 
the condensate population, which in turn drives the system to 
the MI phase at lower interaction strengths and by a strong en- 
hancement of the intrinsic quantum noise signaled in the noise 
correlations. Our analysis demonstrates that these higher or- 
der correlations, which can be measured experimentally in 
time of flight images lillil, provide a unique experimen- 
tal probe for detecting different quantum phases in rotating 
lattices. 

We consider a system of N ultra-cold bosons with mass M 
confined in a uniform ID ring lattice of L sites with lattice 
constant d. The ring is rotated about its axis {z axis) with 
angular velocity Q, . In the rotating frame of the ring the many- 
body Hamiltonian is given by 



H 



2M 



V2 



(1) 

where a is the s-wave scattering length,F(x) the lattice po- 
tential and Lz the angular momentum. Assuming that the 
lattice is deep enough to restrict tunneling between nearest- 
neighbor sites and the band gap is larger than the rotational 
energy, the bosonic field operator, $ can be expanded in Wan- 
nier orbitals confined to the first band $ — "Yli-j O'jWji^)^ 



exp 
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M/j (x). Here Wj (x) 



■ ^ A(x') 

are the Wannier orbitals of the stationary lattice centered at 
the site j, A(x) = 17i x x an effective vector potential and 
dj the bosonic annihilation operator of a particle at site j . 
In terms of these quantities, the many-body Hamiltonian can 
be written, up to onsite diagonal terms which we neglect for 
simplicity, as llsiflHl: 



H = 



e '"d^dj+i + e'"djd]^^ + - 
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Here n 



J — ajCij, 9 is the effective phase twist induced by 

MU^^ and J is 

h ' 



the Gauge field, 9 = /^;+' A(x') • dx' 
the hopping energy; J ^ J dxW* 
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and U the on-site interaction energy: U 

We first carry out a BogoHubov analysis(BA) [15] to un- 
derstand the effects of rotation on the bosonic ring in the 
weakly interacting limit. The BA approximates the field oper- 
ator by a c-number plus small fluctuations a„ = z„ + (5„. 
Replacing it in the many-body Hamiltonian and including 
terms up to second order in 5 yields a quadratic Hamiltonian 
which can be diagonalized by the canonical transformation. 



5n = 



Ulan - vj^a 



For the Hamiltonian described by 



Eq.(l2]i the diagonalization procedure yields: 
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ig and vl = ^e"^'^' 



6fc = ^fc, fc and q integers and 
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-2Jcos[(/)fe -9]+ noU 
2 1 eq + Un 
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with 
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where /i is the chemical potential of the system and 

4Jsin2[%] cos[(^fc ^ e]. no ^ n - -2 



q^o % is the con- 
densate density and n ~ uq the condensate depletion. 

If we write 9 — + ^ with m an integer and < 

A9 < 2tt, the ground state energy corresponds to fc = m 
if < A6' < TT and A: = TO + 1 if TT < A6I < 27r. In 
order words, to lower the energy gained by rotation, the con- 
densate acquires a phase and becomes a current carrying state 
with quasi-momentum (QM) ^k. The quantum number k 
can also be viewed as a measure of the circulation or vortic- 
ity of the gas. Therefore, a change in k from m to m + 1 at 

the critical angular velocity 51^ = ^jv/™2^^2 ^ (i-^- = "')' 
implies the nucleation of a vortex in the system. For n^. 
there is always a unique ground configuration corresponding 
to a macroscopically occupied state with QM equal to ^k. 
On the contrary, at ilc the ground state is not unique as the 
k ~ m and /c = to + 1 solutions are degenerate. At this criti- 
cal rotation, we expect the BA to break down as the underlying 
assumption of the existence of a single macroscopically occu- 
pied mode is no longer valid. Explicit comparison between 
the BA and numerical results, discussed later, shows that the 
BA provides a fair description of the system for small U, and 
away from fie. 

To understand the behavior of the system at critical rota- 
tion, we write the many-body Hamiltonian in terms of QM 



operators, &g = E^Li %c 



-i2-Kqj/L 



,g = 0,.. 



1: 



H=Y^Eqb\hq + ^ J2 W^h,+ 
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J=0 
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q,s,l=0 



Here Eg = —2Jcos[(j)q — 9] are single-particle energies and 
the notation 1 1 1 1 ^ indicates modulo L. In the non-interacting 
limit, at 57c, all the + 1 states given hy \n, N — n) — 



—r^ ^l^^r"^ |0), with n = 0, . . . , iV are degener- 

ate and span the ground state manifold. For small U, we apply 
first order perturbation theory by diagonalizing the Hamilto- 
nian in the degenerate subspace. In view of the absence of 
direct coupling between the different states, the matrix is di- 
agonal and the energy shifts are given by: iji^"* — E.^N + 
^ [{N + n)2 — 3n^ — N]. The states with minimal energy are 
those corresponding to n = and n = A^. To first order in U, 
these two states remain degenerate and higher order perturba- 
tion theory is needed to break the degeneracy. To understand 
the effects of interaction on these degenerate states, one has to 
distinguish between CO and ICO cases. 

In the ICO case, it has been shown[9] that there is no cou- 
pling between | A^, 0) and |0, A^) and hence these states remain 
degenerate for all values of U. Therefore, one can take one of 
these states as the macroscopically occupied mode and use the 
BA to take into account the depletion introduced by quantum 
fluctuations asU/J increases. 

In contrast, in the CO system, there are many different paths 
that couple these two states and their number increases expo- 
nentially with the number of atoms and the sites. The coupling 
leads to the opening of a gap AE. If we write an effective 
Hamiltonian between the two states, we find that 

V^2 = {0,N\H,ff\N,0) = Ai?= (^f (6) 



be- 



ing Hij , the transition matrix elements introduced by the in- 
teraction term of the Hamiltonian and £; the non-interacting 
many-body eigenenergies of the intermediate states. The fac- 
tor — 1) corresponds to the minimum number of collision 
processes necessary to couple the states | A^, 0) and |0, A^) and 
the sum is over all the different paths that generate such cou- 
plings. The ground state of the system becomes : 



ai|0,A^)+a2|A^, 0) 
V2 



(7) 



with 



= —f^ji (see f^). Note that any infinitesimal value 
of U is enough to generate the superposition. As the inter- 
action increases, the states corresponding to other QMs also 
contribute to the ground state and at some point the system 
becomes a MI state. Therefore, at criticality (51 = 51c), the 
ICO and CO respond very differently to the rotation. 

Fig. 1 shows the variation in the spectral gap with inter- 
action for different values of rotation. The critical rotation 
is characterized by the opening of a gap at small U with a 
power-law dependence on U as predicted by Eq. |6l Further- 
more, the compressibility (d^/dh) shows a singular behavior 
at integer fillings(see the inset in Fig.l). Away from critical- 
ity, the system undergoes the usual Mott transition at finite but 
reduced value of U/J. The decrease in the MI threshold can 
be explained by means of a mean field analysis and will be 
discussed later in the paper. 

The complexity inherent to the highly entangled nature 
of the ground state as 51 approaches 51c, and the strik- 
ing different behavior between the CO and ICO system 



3 




0.1 1 10 

U/J 



FIG. 1: (color online)Gap for L = 9 and TV = 9 as a function of 
U/J. The solid (red), dashed (blue), dotted-dashed (green) and dot- 
ted (black) curves are for AS/27r = 0, 0.3, 0.48 and 0.5 respectively. 
The grey line is proportional to (C//J)"'^~^'. The dots correspond 
to Uc calculated from the maximum in the noise correlations ( see 
Fig. 3). The inset shows the variation in the chemical potential (p) 
vs density (n) and [7/ J at f^c. 

can be best illustrated by the first and the second order 
correlations, namely the momentum distribution n{Q) = 
rE™,„e^'"'^^"""^/'^aL«", g = 0,1,...,L - 1 and 
the noise correlations, A(Qi,Q2) = {n{Qi)n{Q2)) — 
{n{Qi)){h{Q2)). 




FIG. 2: (color online) n(Q)/L vs. U/J and A0 for L = 9: Right 
A'^ = 9, Left TV = 8. Black: Q = k (condensate mode ), red: 
Q = fc ± 1 (± depending if ^ vr) 

In Fig. 2 we display the fractional number of atoms in 
the Q — k and Q = k ± 1 modes as a function of Ad and 
U/J, with ± determined by Ad ^ n. In the ICO case there 
is a single macroscopically occupied mode even at ftc, and 
the nature of the ground state is weakly dependent on 6. In 
contrast, for the CO system there is an abrupt redistribution 
of the population of atoms as 6 approaches flc- At criticality, 
for any U/J > 0, the population of the Q = k and Q = 
fc ± 1 modes becomes identical. While for small U, only these 
two modes are occupied, other modes become populated as U 
increases and for U/J ^ 1, irrespective of the 6 value, the 
state of the system becomes an equal superposition of all the 
different QM components. 

The sudden rearrangement of atoms in different QM states 
changes the topology of the n{Q — fc, U, 9) and n{Q — 
k zt 1,U, 6')surfaces: they touch at the critical rotation fre- 
quency illustrating the entanglement. This redistribution of 



momentum, accompanied by an opening of a gap in the spec- 
trum and a non-analyticity in the compressibility ( see fig. 1) 
is reminiscent of the Lifshitz transition in fermionic system 
where a change in the pressure beyond a critical value results 
in an abrupt change in the Fermi surface. In other words, rota- 
tion mimics the effect of pressure and induces a Lifshitz-like 
transition in bosonic atoms confined to a ring shaped geome- 
try. 

The existence of maximal entanglement at the onset to the 
transition can be effectively demonstrated by the noise cor- 
relations which maps the intrinsic quantum noise into an in- 
tensity signal. Away from the critical frequency, the onset 
of the superfluid to MI transition at Uc is signaled by the 
appearance of a maximum in the intensity of A{k,k) ITtIi . 
However, at flc (see Fig. 3), the maximal entanglement oc- 
curs at infinitesimal U, coinciding with the formation of the 
cat state. As a consequence the usual peak at Uc is trans- 
formed to an abrupt jump in the second-order correlations 
at [/ > 0. This jump is coiTectly predicted by Eqs. (Q, 
which yields h{Q) = ^('^'"^Q + ^™+i,q) A(m, m) = 
A(m + 1, m + 1) = — A(m, m + 1) = and zero other- 
wise, in excellent agreement with the numerical simulations. 
Fig. 3 also shows a faster convergence to the asymptotic be- 
havior, A{Q, fc) 26Q^k and n{Q) n, characteristic of a 
system deep in the MI phase with suppressed number fluctua- 
tions lilliTti. 
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FIG. 3: (color online) (a-c) shows A(Q — fc, fc) vs Q — fc and U/ J 
for L = 9 and TV = 9 and d/L = 0, 0.4, 0.5. The plots show that 
A(Q — fc, fc) exhibits a maximum at the onset of the MI for < 9^- 
At 6c a jump at (7 > occurs. In the inset (d) we plot the intensity 
of the maximum intensity reached by A(fc, fc) (triangles right scale) 
and the U/J value at which it takes place (boxes left scale) 

The inset of Fig. 3 illustrates the enhancement of quantum 
correlations as the system approaches ilc- For a given A9, we 
plot the maximum value of A{k, fc) as U/J is varied. Further- 
more, we also show the actual U/ J value at which the max- 
imum is reached. As stated earlier, this value for ^ ilc is, 
related to the onset of the MI transition, and so the plot demon- 
strates a shift in Uc as the the rotation approaches the critical 
frequency. The Gutzwiller ansatz (GA)118] provides further 
insight into the MI transition in the rotating ring lattice. The 
generic GA is based on the assumption that the wave function 
can be approximated as a product of the wave functions at the 
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various lattice sites, I (/3g) = Ilf=o {j2'^=o f^'^ \'^) i) ■ ^^^^e 
\n)i is a Fock states with n atoms at site i. Deep in the MI 
phase, the only relevant states are those characterized by n, 

n + 1 and n—1 atoms per site. Assuming f^^ = ^ i^e"^''-', 

ing the energy E = {Lpa\H\ifiG) with respect to e and e is 
found to vanish at Uc given by: 

Uc = 2Jcoa{27rAe/L){V^+ VnTlf (8) 

where A6 = 2nk/L with fc = to for < < tt and 
fc = TO + 1 for vr < A^? < 2tt. Thus, the GA predicts a 
degenerate ground state at ftc, in agreement with BA. Fur- 
thermore, it also shows a decrease of the critical value of the 
MI transition threshold with increasing rotation by a factor 
of cos{2ttA9/L). The shift is much smaller than that value 
found in our numerical calculations, but this is expected be- 
cause mean field theories only provide a qualitative analysis 
in ID systems. 

For the ICO system, correlations can be computed ana- 
lytically in both small and the large U limit. For small 
U, using Eq. (O and (|4|i, we can determine the n{Q) and 
A((5, Q') . The resulting momentum distribution fi{Q) — 
Lna6k,Q + vf,_Q implies a sharp interference peak at the 
QM of the condensate whose width increases with interac- 
tions due to quantum depletion. On the other hand, noise 
correlations are directly related to the quantity UqVq, the so 
called anomalous average Il6ll which is a measure of the 
many-body scattering matrix and accounts for the modifica- 
tion of binary scattering properties due to the presence of other 
surrounding atoms. The only non-vanishing correlations are 
A(fc, k), A{k, Q), A(Q, Q) and A(2fc - Q, Q) with Q ^ k. 
A{k, Q) ~ ^2u^._qU^._q probes pair excitations of the con- 
densate and has a negative intensity. The anti-correlation can 
be understood as a consequence of the destructive interfer- 
ence between the excited pairs generated when two conden- 
sate atoms collide. A{Q, Q) — A{Q, 2k — Q) = u\_Qvf,_Q 
probe collisions between atoms out-of the condensate; these 
processes are less frequent than condensate collisions and that 
is why they have half of the intensity of A(/c, Q). A{k, k) = 
2 J2q^o '^q'^q^ ^as the larger intensity as its value is twice the 
sum of the intensities of the other A{Q, Q) peaks. Fig. 4 
shows a comparison between analytic results obtained within 
BA approximation and the exact numerical diagonalization 
of Eq. (2) for a finite ICO system. We plot A(to,to), 
A(m + 1, TO + 1) and A(to, to -|- 1) as a function of the ra- 
tio [// J for different rotation velocities. The agreement in the 



small U for a system of only 9-lattice sites is reassuring that 
our analysis based on a finite lattice with few atoms captures 
relevant features. 

In the strongly interacting limit, the behavior of the system 
can be described by mapping the repulsive bosons into a gas 
of ideal fermions 1 19]. For simplicity we will assume n < 1, 
although the mapping to fermions can be generalized to larger 
filling factors 1,20.1 ). In this limit, the single-particle energies 
are given by Eq = — 2Jcos(7r/L(2g — 1) — 9) when N is 
even, and = -2 Jcos(7r/i(2(7)-6') forodd TV [HI. These 
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FIG. 4: (color online) Noise correlation vs. U/J for L = 9 and 
N — 8: Exact (solid lines), BA ( black lines). We use no symbols for 
A(m, m), triangles for A(m, m + l)and boxes for A{m + 1, m + 1). 



equations clearly show that for an ICO system at flc there is 
always one empty energy level with the same energy as the 
Fermi level. Therefore, the two-fold degeneracy of the ICO 
ground state can be seen even in the U oo limit. 

In summary, our study demonstrates the effectiveness of 
experimentally-accessible noise correlations in capturing the 
complexity of highly entangled states, and in heralding a novel 
type of phase transition. Our analysis also shows the appeal 
of noise correlation spectroscopy as a probe of the onset of the 
MI transition in finite systems, where the condensate popula- 
tion decreases gradually and thus the momentum distribution 
does not provide a distinctive signature. Our studies raise im- 
portant open questions regarding the onset and nature of the 
MI at flc as instead of the standard transition from a superfluid 
to a MI here we have a transition from a gapped, fragmented 
and incompressible state to a MI. We hope that our suggestion 
of Lifshitz-like transition in bosonic system will stimulate a 
new line of research in condensed matter community. 
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